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ABSTRACT 

Context. Atomic data is crucial for astrophysical investigations. To understand the formation and evolution of stars, we need to analyse 
their observed spectra. Analysing a spectrum of a star requires information about the properties of atomic lines, such as wavelengths 
and oscillator strengths. However, atomic data of some elements are scarce, particularly in the infrared region, and this paper is part 
of an effort to improve the situation on near-IR atomic data. 

Aims. This paper investigates the spectrum of neutral scandium. Sc i, from laboratory measurements and improves the atomic data of 
Sc I lines in the infrared region covering lines in R, I, J, and K bands. Especially, we focus on measuring oscillator strengths for Sc i 
lines connecting the levels with 4p and 4s configurations. 

Methods. We combined experimental branching fractions with radiative lifetimes from the literature to derive oscillator strengths (f- 
values). Intensity-calibrated spectra with high spectral resolution were recorded with Fourier transform spectrometer from a hollow 
cathode discharge lamp. The spectra were used to derive accurate oscillator strengths and wavelengths for Sc i lines, with emphasis 
on the infrared region. 

Results. This project provides the first set of experimental Sci lines in the near-infrared region for accurate spectral analysis of 
astronomical objects. We derived 63 logigf) values for the lines between 5300A and 24300A. The uncertainties in the/-values vary 
from 5% to 20%. The small uncertainties in our values allow for an increased accuracy in astrophysical abundance determinations. 
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1. Introduction 

Knowing the Milky Way’s chemical composition and its chemi¬ 
cal abundance can lead us to understand how the Milky Way has 
formed and evolved. With time, the amount of heavy elements 
in a galaxy increases as new stars are born and die. Studies of 
elemental composition of stars will give information about the 
environment in which they were born. This, in turn, helps us to 
construct the evolution history of our Galaxy. For this reason, it 
is important to derive accurate abundances of stars. 

Elemental abundances can be derived by comparing an 
observed stellar spectrum with a synthetic spectrum. This 
approach requires atomic data, meaning oscillator strengths 
(/-values) and wavelengths. The measured equivalent width of a 
line is directly proportional to the abundance of the element and 
the/-value. The lack of atomic data makes it difficult to compare 
the observed and synthetic stellar spectra. Thus, an accurate and 
complete set of oscillator strengths and wavelengths are crucial 
for abundance analysis. 

In the periodic table, scandium is between the a-elements, 
e.g. O, Mg, Si, S, Ca, and the iron peak elements, such as 
V, Cr, Mn, Fe. The formation scenarios of a-elements and 
iron peak elements have been studied in detail. Stars that end 
their lives in supernova type II explosions are responsible for 
producing cr-elements, whereas stars that die in supernova type 
la contribute to the formation of iron peak elements and yet 
there is no clear picture of how scandium is formed (|Nissen| 


et al.|20d0} Frangois et al.|2004[|Zhang et al.|2008t[Battistini &] 

Bensby|2015|l. 

Studies of F and G stars may provide insight into the cre¬ 
ation of scandium. However, the abundance varies in these stars. 
There seems to be a disagreement in the scandium abundance 
pattern. In some studies, scandium is overabundant relative to 
iron (|Zhao & Magain|l990t Nissen et al. 2000|l, whereas the 


others found no deviations ( IGratton & Sneden|I99I[ |Prochask'a| 

& McWilliam 2000 1 . Furthermore, scandium is an important 
element for understanding Am and Fm stars, which are over¬ 
abundant in iron peak elements but underabundant in Sc and Ca 
( [Leblanc & AIecian||2008| [Alecian et al.||2013] l. In addition to 
these, some studies cannot be completed owing to the quality of 
oscillator strengths in the infrared region ( [Schmidt et al.|2006l l. 
Improved atomic data with small uncertainties in oscillator 
strengths can help in solving these problems. 

The scandium abundance in the Sun was determined by 
Grevesse et al.|(1996i to be 3.18 + 0.10 dex, and by [Zhang 


et al. P008[ l 3.13 + 0.05 dex and by [Asplund et al. pOOPj l to~be 


3.15+0.04 dex. These values are consistent within the error bars. 
A high quality laboratory measurement of oscillator strengths 
with small uncertainties is a requirement for determining the 
accurate abundance of scandium in the solar photosphere. 

In astrophysical applications, log(gf) values are usually 
used, where g is the statistical weight of the lower level and 
/ is the oscillator strength (/-value). Parkinson et al. (1976i 
derived the experimental log(g'/) values for Sc i lines involving 
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the ground state in the optical region. [Lawler & Dakin| ( j 1989) 1 
measured experimental log(gf/) values for lines in the optical 
region involving the ground state or to the 4s term. For 
other transitions, generally semi-empiric al values of jKuruczj 
( 2009 j l are the favoured ones. In addition, Oztiirk et al. ( 2014| l 
calculated the oscillator strengths of neutral scandium with 
the Cowan code and the quantum defect orbital code (CDOT) 
methods but did not include any near-infrared transitions. To 
our knowledge, there are no experimental oscillator strengths of 
Sc I lines in the infrared region. 

In recent years, infrared spectroscopy has become important 
for observing astronomical objects. Both new generation of 
ground-based telescopes and space telescopes, which work in 
the infrared region, demand complete and accurate atomic data. 
Atomic data of scandium are scarce in the infrared region and 
missing data constrain analysis of a stellar spectrum recorded 
with new-generation advanced instruments and, as a result, limit 
the construction of a galactic evolution model. 

In this work, we measured accurate wavelengths and os¬ 
cillator strengths of Set lines from high-resolution laboratory 
measurements in the infrared region from the upper odd parity 
4p V, 40°, 2 g" and even parity 4s ^P, 2s terms. 

Section 2 describes the branching fraction measurements. In 
addition, this section includes the experimental setup we used 
during our measurements, calibration of wavenumbers, and the 
uncertainty calculations. In section 3 we present our results with 
their uncertainties and a comparison of our results with previous 
studies. 


2. Materials and methods 


2.1. Branching fraction measurements 

The oscillator strength, f-value, is related to the transition prob¬ 
ability for electric dipole transitions by 


/= 1.499 ■ 10 ^'®, ( 1 ) 

gi 

where gf„ is the statistical weight of the upper level, gi the statis¬ 
tical weight of the lower level, A the wavelength of the transition 
in A, and A„/ the transition probability in s * between the upper 
level u and the lower level 1. 

The branching fraction (BF) for a given transition from the 
upper level u to the lower level I is defined as the ratio of the 
transition probability, Aui, of the transition to the sum of all tran¬ 
sitions from the same upper level to all lower states i. The tran¬ 
sition probability is proportional to the transition line intensity; 
therefore. 


BFui = 


Aid 

'Ll Aui 



( 2 ) 


Lifetime measurements are complementary to BF deriva¬ 
tions. The radiative lifetime of the upper level, t,,, is the inverse 
sum of all transition probabilities from the same upper level 
Li = 1/ L'lAui- The transition probability of a line can be de¬ 
rived by combining experimental BFs and radiative lifetimes; 


Aui - 


BFiii 

Tu 


(3) 


We derived transition probabilities, A„/, with the help of (Eq. 3) 
and converted these values to/-values using (Eq. 1) . 

To derive accurate BEs, all transitions from the same upper 
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level should be included. In this work, some upper energy lev¬ 
els had transitions in two different spectral regions recorded with 
different detectors. In such cases, because BFs require only the 
relative intensities, the spectra have to be put on a same relative 
intensity scale by using a normalisation factor. This normalisa¬ 
tion factor was determined from the Sc i lines, which were visible 
in both spectra. 

Eig. 1 shows the partial energy level diagram of Sc i levels. 
In this figure we have marked the transitions that appeared in 
our spectra. Transition lines from the same upper level were pre¬ 
dicted from the IKuruczj p009|l database and from [Ben Ahmed) 


& Verges (1977 i. These lines were identified in our spectra by 
using FTS analysis software GEit ( [Engstrom[1998)[20T4| l. 

Sc I has only one stable isotope, "^^Sc, with a nuclear spin 
1= 7/2, so Set lines show hyperfine splitting (hfs). Interaction 
between the nuclear spin and the total angular momentum splits 
energy levels into the hyperfine structure. This splitting appears 
as a resolved line splitting or a line broadening in a spectrum. 
Unsuccessful treatment of hfs leads to incorrect abundances, as 
much as several tenths of a dex ( [Nielsen et al.|20d0| . For lines 
with negligible hfs, one can fit a Voigt profile to the line and 
measure the central wavenumber and intensity. However, this 
is not the case for hfs. Having hfs gives rise to a difficulty in 
fitting a line profile since every fine-structure transition appears 
as a broadened or multi-component feature. For this reason, we 
defined the position of lines by the centre of gravity and instead 
of fitting a Voigt profile, we measured the integrated intensity 
under each measured line with the help of GFit. We applied this 
method for the cases in which we have a single line profile i.e. 
to test the accuracy. In all of these cases, fitting a line profile and 
taking an integrated intensity gave the same result. 


2.2. Experimentai setup 

A water-cooled hollow cathode discharge lamp (HCL) was used 
as a source to produce scandium atoms. The HCL consists of 
a glass tube, anodes on each side, and an iron cathode in the 
middle. A small solid sample of scandium was placed in the 
iron cathode. The diameter of the inner cathode was 7 mm and 
the distance between anode and cathode was 20 mm. The light 
source was run with a pressure of 1.0 Torr (1.3 mbar), with 
different applied currents, ranges from 0.10 A to 0.60 A, and 
with argon or neon as a carrier gas. The best condition, i.e. the 
strongest lines, for the measurements was achieved with argon 
gas and with applied current of 0.60 A. The argon lines were 
used as reference lines for wavenumber calibration as described 
in section 2.3. 

The Set spectra were recorded using the high resolution 
Fourier transform spectrometer (FTS), Bruker IFS 125 HR, at 
the Lund Observatory (Edlen Laboratory). This instrument has 
the wavenumber range of 50000 - 2000 cm“' (200 - 5000 nm) 
and the maximum resolving power of 10^ at 2000 cm '. The 
resolution of the FTS was set to 0.02 cm ' during measure¬ 
ments. We recorded the spectra with the indium antimonide 
(InSb) and silicon (Si) detectors, which are sensitive to different 
spectral regions. 

The optical elements in FTS contribute to the response 
function of the instrument. Observed line intensities differ from 
the intrinsic line intensities, due to the wavelength dependent 
transmission of the spectrometer and optical elements. It was 
therefore required to find the response functions of two different 
detectors. After determining the response function, one can 
acquire the intrinsic line intensities by dividing the observed 
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Fig. 1: Partial energy level diagram of Set levels showing the observed transitions in the measurements. 


line intensities with the response function. A common way to 
determine the response function is to use an intensity calibrated 
reference. We used a tungsten filament lamp to calibrate the 
intensities of Sc i lines. This lamp in turn has been calibrated by 
the Swedish National Laboratory (SP) for spectral radiance in 
the range 40000 - 4000 cm“' (250 - 2500 nm). If the radiance 
of the reference is known, one can determine the response 
function. Fig. 2 shows the response function of the instrument 
with Si and with InSb detectors. 

We recorded the spectrum of the tungsten lamp immediately 
before and after every scandium recording to verify the response 
of the instrument did not change during the measurements. For 
this reason the HCL and tungsten lamp were placed at the same 
distance from the FTS and a rotating mirror was used in order 
to change the light source (HCL or the Tungsten lamp) without 
moving the lamps. There were no changes in the spectra of the 
tungsten lamp taken before and after the HCL measurements, 
therefore one response function for each detector was used 
for calibration. Using intensity calibrated spectra, relative 
line intensities at the different regions were connected with a 
normalisation factor. This was done by using the Sc i lines in the 


overlapping region of the two detectors, see Fig. 2. 

The spectra of Sci lines were recorded with different 
currents. Higher current values help weaker lines to be visible, 
but introduce self-absorption. Self-absorption has a negative 
effect on line intensities e.g. it changes the observed line profile 
and gives an incorrect intensity. This in turn affects the BF mea¬ 
surements, see Eq. 0. The effect of self-absorption depends on 
the number density of atoms in the plasma and intensity of the 
line. The number density of the plasma increases by increasing 
the current applied to the HCL. The higher concentration of 
scandium ions in the plasma results in more lines affected by 
self-absorption. We performed a set of measurements with 
varying currents to determine if Sci line intensities were 
affected by self-absorption at higher currents. Then we plotted 
normalised integrated intensity ratios of lines from the same 
upper level versus applied current. Integrated intensity ratios 
between lines from the same upper level must be same for 
different applied currents if there is no self-absorption. In our 
measurements some lines showed effects of self-absorption. In 
these cases, we extrapolated integrated intensity ratios to zero 
current, where self-absorption is assumed to be zero, and used 
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Fig. 2: Instrument response function with InSb (solid line) and 
Si (dashed line) detectors, respectively. Each curve is normalised 
to its highest point. 


these values to obtain correct relative intensities. ISikstrom et al.l 
( 2002| l tested this method by comparing the self absorbed Fe 
I lines intensity ratios at zero current with the ratios derived 
by Blackwell et al. \\919\ using absorption measurement. 


They found that the extrapolated intensity ratios agree with 
the Blackwell et al. ( 1979| l intensity ratios. The correction of 
self-absorption introduce uncertainties in the BF determination. 


2.3. Wavenumber calibration 


The HeNe laser in the FTS determines the displacement of the 
FTS moving mirror and gives a multiplicative wavenumber scale 
to the spectrum. However, a non-parallel alignment with the in¬ 
coming light beam adds a shift to the measured wavenumber. In 
addition to this difference, having a finite size aperture may lead 
to wavenumber shifts ( |Feamer & Thorne|1988| l. The shift intro¬ 
duced by these effects is multiplicative. With a correction factor 
this shift can be calibrated; 


CTfrue ~ rr^pp^y^fifit -H k^yy). (4) 

In Eq. 0, (Ttrue denotes the corrected wavenumber, (Tapparent is 
the measured wavenumber, and k^yy is the correction factor. We 
used the wavenumbers of Ar i lines from Whaling et al. (2002 1 , 
including the multiplicative correction by Sansonetti (20071, as 
a wavenumber standard for the wavenumber calibration of Sc i 
lines. 

In principle, using only one calibration line is enough to fix 
the multiplicative constant, but several calibration lines can be 
employed for the purpose of increasing the wavenumber accu¬ 
racy. We calculated a correction factor for each Ar i line and used 
the weighted mean of these factors. 

We combined the errors in the determination of the multi¬ 
plicative factor and the errors from centroiding the position of 
each line to estimate an absolute uncertainty of the observed 
wavenumbers of M(cr)=+0.001 cm"^ 


2.4. Uncertainties 

The contribution to the f-value uncertainties comes from the 
uncertainty of the radiative lifetime and the uncertainty of the 
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branching fractions. The branching fraction uncertainties con¬ 
tain the uncertainty of the intensity calibration lamp and the un¬ 
certainty that arises from using two different detectors, i.e. un¬ 
certainty of the normalisation factor and the contribution from 
the self-absorption correction. 

The uncertainty of the line intensity (or the integrated in¬ 
tensity) and the uncertainty of the radiative lifetime are uncor¬ 
related. Although BFs are dependent on each other, we treated 
them as independent in the uncertainty calculations. By includ¬ 
ing all the uncertainties from different effects described above, 
one can determine the total uncertainty of the BFs as specified 
in |Sikstr6m et al.|P002| l, 


(^) 


= (1 - (BF),f 



+ (BF)] 

i*k(inP) 

+ (BF)] 

ji^KinQ) 




(5) 


In the first term on the right-hand side of Eq. 7, denotes 
the branching fraction of the line in question in the spectral re¬ 
gion of the detector P and m( 4) denotes the uncertainty in the 
measured intensity of the line k. The second term includes the 
BFs and measured intensity uncertainties of the other lines from 
the same upper level in the spectral region of detector P and 
the uncertainty of the calibration lamp u{cj). The last term con¬ 
tains the BFs and measured intensity uncertainties of the lines in 
the spectral region of detector Q, the uncertainty of the calibra¬ 
tion lamp, and the uncertainty of the normalisation factor, u{nf), 
between the spectral regions. The uncertainty of the calibration 
lamp, u{cj), is 7% and the uncertainty of the normalisation factor 
is u{nf) = 5%. We used the uncertainties of the integrated inten¬ 
sity, u(I), from GFit and the relative uncertainty varies between 
0.002% for the strong lines and 10% for the two weak lines and 
self-absorbed lines. The cases in which we had self absorbed 
lines, we included an uncertainty from self absorption correc¬ 
tion as well. This uncertainty varies between 1% and 9%. From 
Eq. 3, the uncertainty of the/-values or the transition probability 
becomes 

ft j I A, j i BE j 



where m(t) is the uncertainty of the radiative lifetime of the up¬ 
per level. The derived uncertainties of/-values range between 
5% for the strong lines and 20% for the weak lines or for the 
lines with high uncertainty in their radiative lifetime. 

Two different sets of radiative lifetimes were used, one set of 
experimental lifetimes and one set of semi-empirical lifetimes. 


The uncertainty of the experimental lifetimes is 5% (Marsden 
|et al.|19^ . For semi-empirical values, we compared the semi- 
empirical lifetimes ( Kurucz|2009| l of the levels that have experi¬ 
mental lifetime values and examined how they differ from the ex¬ 
perimental lifetime values. Based on this comparison we adopted 
20% relative uncertainty for the semi-empirical lifetimes. 


3. Results and conclusion 

We calculated the BFs from observed line intensities from Eq. 
2. All the lines from the same upper level should be included to 
have accurate BFs. For this purpose we recorded several spectra 
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Fig. 3: Comparison between log(g'/) values of this work and the 
semi-empirical |Kurucz ( |2009| l log(gf /)^09 values. Upper panel 
(a) includes the lines with semi-empirical lifetimes, whereas the 
lower panel (b) includes lines with experimental lifetimes. 


with different detectors, which are sensitive to different regions. 
We used the lines in the overlapping regions of two different de¬ 
tectors to put lines on the same relative intensity scale. Some of 
the weak lines were not observed in the spectra. In these cases we 
used theoretical transition probabilities to estimate the residual 
value, i.e. the missing BF from the unobserved lines, for these 
levels. The values are less than one per cent for all levels, except 
from the 4s ^Si /2 level. Derived log(g/) values for the lines to¬ 
gether with their uncertainties and the residual values of missing 
lines are given in Table 1 and Table 2. 

In addition to using different detectors, the HCL was run 
with different discharge currents. This helped us to find out if 
the strong lines were affected by self-absorption. In the cases in 
which there were lines affected by self-absorption, we corrected 
the measured intensities of these lines. This in turn added an ex¬ 
tra uncertainty to the BF of the affected line. 

We derived transition probabilities, Aui, by combining BFs 
and radiative lifetimes. Where available the radiative lifetimes 


Fig. 4: Comparison between \Qg(gf) values of this work de¬ 
rived from the experimental radiative lifetimes and the previous 
experimental \ogigf)pRK values of Parkinson et al.|p976| l and 
\og{gf)Lw values of [Lawler & Dakin U989| l, in subplot (a) and 
(b), respectively. These studies do not include infrared transi¬ 
tions. 


are taken from the measurements of |Marsden et al.| ( |1988| l and 
the rest of the lifetimes are semi-empirical values from |Ku-| 


rucz 


(2009|. The dominating factor in the uncertainty of oscil¬ 
lator strengths comes from the radiative lifetimes, resulting in 
higher uncertainties for the oscillator strengths derived from the 
semi-empirical lifetimes. They could be improved with labora¬ 
tory measurements. 

We derived oscillator strengths for 63 lines from 4p 


4t:?o 4r^(> 2r^D 2r)0 Ar'o 

^ -i 12 , 512 , 112 , 912 y ^ 1/2,3/2,5/2,7/2, ^ 3/2,5/2, ^ 3/2, ^ 5/2,7/2,9/2, 

^G° 9 / 2 , and 4s ^Pi/ 2 , 3 / 2 , ^Si /2 levels with uncertainties in the/- 
values between 5% for the strong lines and 20% for the weak 
lines or for the lines with high uncertainty in their radiative life¬ 
time. These results together with their branching fractions, tran¬ 
sition probabilities, and previously published log(gf/) values are 
presented in Table 1 and Table 2. The first column of the previ¬ 
ously published \Qg{gf) values is the semi-empirical calculations 
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of [Kurucz (2009 1 and the other two columns are experimental 
values of[ ^rkinson et al.| ( |1976| l and [Lawler & Dakin| ( |1989| l. 
These previous measurements only included lines for the optical 
region and to our knowledge this is the first time that the exper¬ 
imental logigf) values of Sc i lines in the infrared region were 
measured. Table 1 shows our results derived from the experi¬ 
mental radiative lifetimes and Table 2 from the semi-empirical 
radiative lifetimes. 

In Figs. 3 and 4, we graphically compare our data with pre¬ 
viously published studies. Fig. 3 shows our results compared 
with Kurucz (20091 semi-empirical \og{gf)Kog values and Fig. 
4 presents our results compared to experimental \og(gf)pRK of 
[Parkinson et al.| ( ri976[ ) and log(g/)LW' [Lawler & Dakin] ( [1989| l. 
The \og{gf) values in Fig. 3a were derived from the semi- 
empirical radiative lifetimes, and in Fig. 3b they were derived 
from the experimental values. As seen in Figs. 3a and 3b, the 
\og{gf) values measured in this work agree with the semi- 
empirical values for high logigf) values. For lower values the 
scatter is larger. This can be explained by the difficulties in cal¬ 
culating spin-forbidden lines theoretically, uncertainties in the 
semi-empirical calculations and by the fact that weaker lines 
have larger uncertainties than strong lines in the experiments. 
In Fig. [TO and Fig. 4a we compared our results derived 


from the experimental lifetimes with the previous experimen¬ 
tal \ogigf)pRK values of Parkinson et al.|p976| l and \og{gf)Lw 
values of Lawler & Dakin ( 1989| l in the optical region. Com¬ 
parisons with experimental values show very good agreement, 
especially with [Lawler & Dakin[ ( [1989[ ). This supports that our 
results are accurate in the infrared region. The small uncertain¬ 
ties for the oscillator strengths presented in this work allow for 
improved stellar abundances. 
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